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Abstract 

The like sign dileptons provide the most promising signature for superparticle search in 
a large category of i?-parity violating SUSY models. We estimate the like sign dilepton 
signals at the Tevatron collider, predicted by these models, over a wide region of the MSSM 
parameter space. One expects an unambiguous signal upto a gluino mass of 200 — 300 GeV 
(> 500 GeV) with the present (proposed) accumulated luminosity of ~ 0.1 (1) fb _1 . 
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Introduction 



The canonical signature for superparticle search at hadron colliders is the missing-p^ 
signature, which follows from i?-parity conservation [1]. The latter implies that the super- 
particles are produced in pair; and the pair of lightest superparticles (LSP) resulting from 
their decays are stable. In the SUSY models of current interest the LSP happens to be 
the lightest neutralino Z\. Being a weakly interacting particle, the LSP escapes detection 
like the neutrino. The resulting imbalance of visible transverse-momentum constitutes the 
missing-pr signature for superparticle production. However, there is a growing interest now 
in possible signatures for superparticle production in i?-violating SUSY models, since there 
is no compelling reason for .R-parity conservation in supersymmetry. It is usually invoked to 
ensure proton stability via the conservation of lepton and baryon numbers, L and B. How- 
ever proton stability requires the conservation of L or B, but not necessarily both. Thus one 
can have two types of i?-violating SUSY models, corresponding to B or L violation, which 
are consistent with proton stability. The former implies LSP decay into a baryonic channel, 
which has a large QCD background. On the other hand the latter implies LSP decay into a 
leptonic channel, which can serve as a viable signature for superparticle search. In particular 
the Majorana nature of the LSP (Zi) implies that the LSP pair decay into like and unlike 
sign dileptons with equal probability. This leads to a like sign dilepton (LSD) signature for 
superparticle production, which has very little background from the standard model (SM). 
The dilepton data from an early Tevatron run [2], corresponding to an 
integrated luminosity of ~ 4 pb" 1 , was analysed in [3] in terms of the above signature. 
This led to a lower bound on gluino and squark masses 

m~ g ,mq > 100 GeV, (1) 

in the /^-violating SUSY model. At present each Tevatron experiment has accumulated data, 
corresponding to an integrated luminosity of ~ 100 pb _1 . Moreover this is expected to go up 
by another order of magnitude following the main injector run (run II). Thus it is imperative 
to explore the prospect of extending the search to higher superparticle masses in the above 
mentioned i?-violating SUSY model using these data. The present work is devoted to this 
exercise. 

We shall work within the framework of the minimal supersymmetric standard model 
(MSSM), which is characterised by relatively few parameters [1]. For extending the analysis 
to rrig > 100 GeV, one has to include two additional features, which were not considered in 
[3]. Firstly one has to consider the cascade decay of gluino via the chargino and the heavier 
neutralino states, which dominate over its direct decay into the LSP for m g > 100 GeV. 
Consequently the result will depend not only on the gluino mass, but on the two other 
MSSM parameters as well - i.e. the higgsino mass parameter (/i) and the ratio of the 
two higgs vacuum expectation values (tan/3). Secondly one has to include the contribution 
from the electroweak production of these chargino (^1,2) and neutralino (-Zi_ 4 ) states [4]. 
As we shall see below, this contribution dominates over the gluino decay contribution for 
m- g > 200 GeV. 

In the next section we summarise the i?-violating decays of LSP, which will be relevant 
for our analysis. In the following section we shall briefly discuss the processes contributing 
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to LSP production, i.e. the gluino cascade decay and the electroweak production of the 
chargino and neutralino states. The resulting signal cross-sections will be presented in the 
next section. We shall conclude with a summary of our results. 

LSP Decay in i?-violating SUSY Models 

There can be explicit [5,6] as well as spontaneous [7] breaking of i?-parity; but only the 
former is phenomenologically viable in the MSSM. Therefore we shall concentrate on explicit 
-R-parity breaking, where the LSP decay arises from one of the following Yukawa interaction 
terms in the Lagrangian: 

Cft = \ijddjGk + X' ijk £iqjd k + \" jk didjU k (2) 

plus analogous terms from the permutation of the supertwiddle. Here £ and e (q and it, d) 
denote the left-handed lepton doublet and antilepton singlet (quark doublet and antiquark 
singlet) and i,j, k are the generation indices. Evidently the first two terms violate L and 
the third one violates B conservation. The A and A" couplings are antisymmetric in the 
first two indices, so that there are 9 independent ones of each. Together with the 27 A' 
couplings, there are 45 independent Yukawa coupling terms. In analogy with the standard 
higgs Yukawa couplings, one expects a hierarchical structure for these additional Yukawa 
couplings as well [5,6]. The decay channel of LSP is determined by the leading Yukawa 
coupling term. 

The leptonic decay channels of LSP correspond to one of the A or A' couplings being the 



leading one, as shown below. 

Kjk ■ Zi -»• Tqq'ivrqq'), (3) 

Ksk,23k ■ Zi -»■ "eW, (4) 

K jk (h3^3) : Z^iqq'^qq'), (5) 

Ai33,233 : Z x -> £u T f(veTf), (6) 

A 123 : Z x ^lv t ,f{y t f!f), (7) 

A 3 11,322,312,321 = Z x -> TV($ \v T if) , (8) 

A m ,i22 : Z x ^lv v ?{y t W), (9) 



where i denotes e, /x and each of the above final states represents the corresponding charge 
conjugate state as well. Note that in 15 of these 36 cases, represented by (3) and (4), there 
is no e or /i in the final state. The £tq' final state in (4) is kinematically inaccessible for 
the Z\ mass range of interest. Thus the leptonic decay channels of our interest correspond 
to the remaining 21 cases, represented by eqs. (5-9). The corresponding squared matrix 
elements are given in [8,9]. For the A couplings, the pairs of final states shown in eqs. (6- 
9) have a branching fraction of 50% each, assuming a common selectron mass. For the A' 
couplings of eq. (5), however, the branching fraction for the Z\ — > £qq' decay is sensitive to 
the composition of Z\. This is shown in Table I for different values of the MSSM parameters, 
where we have assumed a common sfermion mass = rriq ^> rag . 

It may be noted here that 15 of these 21 cases, represented by eqs. (5-7), lead to no more 
than two leptons in the decay of the LSP 
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pair. Moreover for the 4 cases represented by eq. (8) the dilepton final state dominates 
over the 3-4 lepton state, as we shall see below. Thus in 19 of the 21 cases of interest, the 
LSD channel is the most viable channel for superparticle search. Only for the last two cases, 
represented by eq. (9), the trilepton final state dominates over the LSD [10]. Evidently they 
represent the most favourable case for the multilepton signal of superparticle production. 
We shall not present the signal cross-section for this case, since it has been covered in [10]. 

We shall conservatively assume the leading /^-violating Yukawa coupling to be signifi- 
cantly less than 1, so that the pair production of superparticles and their decays into LSP 
are not affected. In this case the signal does not depend on the value of the Yukawa coupling 
as long as it is large enough for LSP decay inside the detector, i.e. ^ 10~ 5 [3,6]. 

Production and Decay of Superparticles into LSP 

The MSSM implies rriq ^ m g . In estimating the gluino and chargino/neutralino cross- 
sections we shall conservatively assume rriq to be significantly larger than m g . In that case 
the cross-sections are insensitive to rriq [11]. The results presented below are obtained with 

rriq = 2m g . (10) 

The dominant processes for gluino production are the leading order QCD processes [12] 

gg(qq)^gg. (11) 

In order to discuss the cascade decay of gluino, a brief summary of the chargino/neutralino 
sector is in order. The masses of the SU(2) and U(l) gauginos, M 2 and M 1: are related to 
the gluino mass in the MSSM [1], i.e. 

M 2 = m- g ~ 0.3m- g , (12) 

M x = \ tan 2 9 W M 2 ~ 0.5M 2 . (13) 

The physical neutralino states -Zi-4 are mixtures of these two gauginos and the two neutral 
higgsinos. Similarly the physical chargino states W^ 2 are mixtures of the charged SU(2) 
gaugino and the charged higgsino. Their masses and compositions are obtained by diagonal- 
ising the corresponding mass matrices [1,13]. They are functions of m g , \i and tan/3 [9,14]. 
It is important to include the QCD correction factor, which relates this running gluino mass 
with its physical (pole) mass, i.e. [15] 



m g (pole) = rrig(mg) 



i + 42a " 
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(14) 



The signal cross-sections are presented below in terms of this pole mass. 

It is worth mentioning here that for most of the MSSM parameter space of our interest 
the higgsino mass parameter ii is > Mi,M 2 . Consequently the lightest neutralino (Zi) is 
dominated by the U{\) gaugino component, while the second lightest neutralino [Z 2 ] and 
the lighter chargino {W\) are dominated by the SU{2) gaugino. Thus their masses roughly 
correspond to 

mz, - M u (15) 
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m z 2 - m m - M 2- ( 16 ) 
Moreover, the gluino decays dominantly into these states, i.e. 

g ™> gg'W,, (17) 
g qqZ 2 , (18) 

g qqZ u (19) 

where the larger branching fractions into W\ and Z 2 reflect the larger SU (2) gauge coupling 
relative to the U (1) [9,14]. Of course our results are obtained with exact values of the masses 
and branching fractions, which show significant deviations from the approximate formulae 
(15-19) over parts of the parameter space. 

In addition to the above, one has to consider the electroweak processes for chargino/neutralino 
production [4,16] 

qq -> W^WiZ^WfWi . (20) 

These are dominated by the s-channel W and Z exchanges. In spite of being electroweak 

processes they dominate over the QCD process (11) for m~ g > 200 GeV because of the 
relatively low W\ and Z 12 masses. 

Finally the W\ and Z 2 coming from (17), (18) and (20) decay into the LSP (Z ± ) via W 
and Z exchanges, i.e. 

W, qq'Z^Zj, (21) 

Z 2 qqZ x {trZ x ). (22) 

It may be noted that the leptonic decays of W\ and Z 2 have branching fractions of 0.22 and 
0.06 respectively. The former can give rise to a LSD signal for gluino production via (11) 
and (17) in the i?-conserving SUSY model. Indeed this signal is expected to be as good as 
the canonical missing-p T signal at the LHC energy [9,14]. For the Tevatron energy, however, 
the size of this LSD signal is rather small [17]. On the other hand, one expects a significant 
contribution to the LSD signal in the i?-violating SUSY model from the cross-term, where 
one of the leptons comes from the W\ (or Z 2 ) decay. This contribution is included in our 
estimate of the signal cross-section. 

The cross-sections presented below are calculated for the Tevatron collider energy of 

= 1.8 TeV, (23) 

using the MRSD'_ structure functions [18], with the QCD scale chosen as the sum of the 
produced superparticle masses. 

Results and Discussion 

We have calculated the signal cross-sections as functions of gluino mass for 

tan = 2, 10 and /i = -100, -200, -300, +300 GeV. (24) 
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A large part of the interval —100 < \x < 300 GeV is excluded by the LEP data for the 
low gluino mass range of our interest [19], while one expects no significant change in the 
cross-sections beyond \fi\ = 300 GeV. 

Fig. 1 (a,b,c,d) shows the cross-sections for gg, W\Z 2) W^W^ and W\Z\ production at 
tan (3 — 2 and the four values of /i mentioned above. One can clearly see the electro-weak 
processes W\Z 2 and W{W\ overtaking the QCD process gg for m g > 200 GeV at negative 
values of fi (Fig. la,b,c). The rapid increase in the W±Z 1 cross-section as m g goes down to 
150 GeV reflects production via on-shell W as +ttl^ i falls below myy. This enhancement 
is visible for all the three electroweak processes at ji = 300 GeV (Fig. Id), as the Zi, Z 2 and 
W\ masses are all small for positive /i. This is why the region m g < 250 GeV is excluded 
by LEP data for positive /i. It should be noted, however, that the electroweak processes 
dominate over the QCD process of superparticle production for any gluino mass at positive 

Fig. 2 (a,b,c,d) shows the cross-sections for the three electroweak processes at tan/? = 10. 
They are significantly larger than the previous case at /i = — 100 GeV (Fig. 2a). This is due 
to a drop in the Zi, Z 2 and W\ masses in going from tan/5 = 2 to 10 at negative /x, which is 
most significant at fi = —100 GeV. Note that all the three cross-sections in Fig. 2a shoot 
up for m g ~ 200 GeV due to on-shell W and Z exchanges as the corresponding thresholds 
fall below the W and Z masses. For the same reason, however, this mass range is ruled out 
by the LEP data. On the other hand the cross-sections are clearly smaller than the previous 
case at /i = +300 GeV (Fig. 2d). This reflects an increase in the above masses in going 
from tan/9 = 2 to 10 at positive fi. Finally the suppression of W±Z 1 relative to the W\Z 2 
and WiWi cross-sections at fi = —200, —300 and +300 GeV is due to the decoupling of the 
SU(2) gaugino component from Z\. All these features can be checked with the masses and 
compositions of these particles listed e.g. in [9]. 

Fortunately the wide variation of the chargino/neutralino cross-sections with /z and tan/3 
parameters does not reflect in the resulting LSD signals. Lower chargino/neutralino masses 
correspond to softer decay leptons; and the resulting reduction in the detection efficiency 
compensates for the rise of the corresponding cross- sect ion. Moreover the LSP (Zi) coming 
out from on-shell W and Z decays carry very little p?, so that its decay lepton seldom passes 
the required px cut. Thus the sharp peaks seen in some of these cross-sections at low m g 
has little effect on the resulting LSD signals, as we see below. 

The LSD signal coming from the above superparticle decays has been estimated with the 
following p T rapidity and isolation cuts on each lepton [2] : 

p e T > 15 GeV, \rj e \ < 1, E T C < 5 GeV, (25) 

where the last quantity refers to the transverse energy accompanying the lepton within a 
cone of AR = ^Ar/ 2 + A</> 2 = 0.4. Our estimates are based on a parton level Monte- 
Carlo programme, which should be adequate for the leptonic signal. The only place where 
hadronisation could play a significant role is in the isolation cut. We have checked that the 
efficiency factors obtained by our parton level programme for the isolation as well as the pt 
and rapidity cuts agree with to those of IS A JET calculation [20] . 

The LSD background from the standard model has been calculated in [10]. The total 
background cross-section is only 2.4 fb, coming from WZ (2.1 fb) and ti (0.3 fb) [21]. This 
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corresponds to a quarter of an event for the accumulated luminosity of ~ 0.1 fb -1 in the 
present Tevatron run, going upto 2.4 events for the projected luminosity of 1 fb -1 with the 
main injector. Thus the discovery pottential of the LSD signature is expected to be limited 
by the signal size rather than the SM background [10]. 

Fig. 3(a,b,c,d) shows the LSD signal cross-sections resulting from gg,WiZ 2 ,WiW{ 
and W\Z X production at tan/? = 2 and \x = —100,-200,-300 and +300 GeV, assuming 
the A' coupling of eq. (5) to be the leading /^-violating Yukawa coupling. The relevant 
branching fractions for Z\ — > iqq' decay are given in Table I. We see from Fig. 3 that the 
LSD signal cross-sections are less sensitive to ji than the corresponding raw cross-sections 
of Fig. 1 as remarked earlier. The WiZ 2 and H^Wf contributions dominate over gg at 
m g > 200 GeV through out the negative \x region, while the W 1 Z 1 contribution is very 
small. At /i — +300 GeV, the electroweak contributions dominate over gg for all values 
of nig. Note that the size of the signal cross-section here is quite similar to that in the 
negative \i region. This is somewhat accidental as the larger cross-section at fj, — +300 GeV 
is compensated by a smaller branching fraction for the Z\ — > iqq' decay (Table I). 

A brief comment on the effect of kimatic cuts is in order. At m g = 300 GeV and negative 
\x the suppression factors from the lepton p T , rapidity and isolation cuts are around 5, 2 and 
2 respectively, resulting in an overall suppression factor of ~ 20. This goes up (down) by a 
factor of 2 at m g = 200 (400) GeV. The corresponding suppression factors at fj, — +300 GeV 
are about twice as large because of the low chargino/neutralino masses. 

Fig. 4 shows the net LSD signal cross-sections for different choices of the leading R- 
violating Yukawa coupling at tan/5 = 2. Let us first consider the negative \x region (Fig. 
4a,b,c). For the A' and the unfavourable A couplings of eqs. (5) and (6) one expects a LSD 
signal cross-section > 100 fb for m g = 200 GeV. For the more favourable A couplings of eqs. 
(7) and (8) the LSD cross-section remains > 100 fb upto m g = 300 GeV. This corresponds 
to at least 10 isolated LSD events for the current Tevatron luminosity of 0.1 fb -1 . Thus the 
current CDF data is capable of probing for /^-violating SUSY signal upto a gluino mass of 
atleast 200 GeV in the former case and 300 GeV in the latter. With the expected luminosity 
of ~ 1 fb" 1 at the main injector run, the probe can be extended upto 500 GeV in the former 
case and ~ 600 GeV in the latter. Turning to the jj, = +300 GeV region (Fig. 4d), one sees 
that the A' coupling predicts no viable LSD signal for the current luminosity of 0.1 fb -1 . 
However the A coupling of eq. (6) predicts a viable signal for m g = 300 GeV, while for 
the more favourable couplings of eqs. (7) and (8) it remains viable upto 500 GeV. With 
a luminosity of 1 fb -1 of course one expects a viable signal upto 500 GeV even for the A' 
coupling, while for the A couplings the discovery limit is significantly larger. It should be 
noted here that the four A couplings of eq. (8) can lead to 3 and 4 lepton final states as well. 
For this case the trilepton signal, corresponding to 3 leptons passing the kinematic cut (25), 
is also shown in Fig. 4. It is comparable to the corresponding LSD signal in the negative \x 
region. But it is relatively small at /i — +300 GeV, as the low chargino/neutralino masses 
make it harder for the 3rd lepton to pass the kinematic cut. The 4 lepton signal (not shown) 
is negligible. 

Finally Fig. 5 shows the signal cross-sections for different choices of the leading R- 
violating Yukawa coupling at tan/3 = 10. The most noticable feature in this case is the drop 
in the signal cross-section for the A' coupling ad, fx — —100 GeV (Fig. 5a). This is due to the 
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fall in the Z\ — > £qq' branching fraction shown in Table I [22]. Apart from this the signal 
cross-sections are generally insensitive to the choice of \x. With the current luminosity of 
~ 0.1 fb _1 , one expects a viable LSD signal for m g = 200 GeV for the A' and the unfavourable 
A couplings of eqs. (5) and (6), while for the more favourable A couplings of eqs. (7) and (8) 
it remains viable upto m g = 400 GeV. With a luminosity of ~ 1 fb -1 , the discovery limit 
goes upto m g = 500 GeV in the former cases and much beyond in the latter. Note that for 
the A couplings of eq. (8) one also expects a trilepton signal of similar size. 

It may be added here that for the most favourable A couplings of eq. (9), the multilepton 
signals studied in [10] are somewhat larger than the ones shown above. Although these 
signals were studied in [10] only for one set of ji and tan/3, they are expected to remain 
comfortably large for other values of these parameters as well. 

Summary 

The multilepton signals provide the most viable signature for a large category of i?-parity 
violating SUSY models, where the LSP undergoes leptonic decay via one of the Yukawa 
couplings of eqs. (5-9). In all but 2 of these 21 cases the like sign dileptons constitute the 
dominant SUSY signal. We have done a systematic analysis of these LSD signals at the 
Tevatron collider energy covering a wide range of the MSSM parameters. The contributions 
from the gluino cascade decay as well as the electroweak production of chargino/neutralino 
pairs are taken into account. For the A' and the unfavourable A couplings of eqs. (5) and 
(6) one expects an unambiguous LSD signal upto a gluino mass of at least 200 GeV with 
the current accumulated luminosity of ~ 0.1 fb -1 . For the projected luminosity of ~ 1 fb -1 
from the main injector run, the signal is expected to remain viable upto a gluino mass of 
500 GeV. For the more favourable A couplings of eqs. (7) and (8) one expects viable signals 
upto a gluino mass of 300 GeV and at least 600 GeV for integrated luminosities of 0.1 and 
1 fb -1 respectively. It may be added here that the most favourable A couplings of eq. (9), 
for which the trilepton signal dominates over the LSD, has been already investigated in ref 
[10]. The signal cross-section in this case is somewhat larger than above. 

At present only the CDF experiment can probe for the LSD signal. But the D0 experiment 
can also probe for this signal in the main injector run, since it is scheduled to install a central 
magnet for lepton charge identification. Finally the increase of the energy from 1.8 to 2 TeV 
in the main injector run will push up the corresponding discovery limits somewhat higher. 
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Table I 



The Branching Fraction for LSP decay [Z\ — > £qq') via the R- violating Yukawa coupling 
Kjk 3 ) 



rriq 


n 


tan = 2 


tan/? = 10 


tan (3 = 30 


250 


-100 


.87 


.54 


.21 




-200 


.85 


.63 


.46 




-300 


.82 


.61 


.50 




+300 


.15 


.24 


.36 


200 


-100 


.85 


.40 


.18 




-200 


.82 


.57 


.43 




-300 


.77 


.58 


.49 




+300 


.126 


.28 


.38 


300 


-100 


.81 


.27 


.16 




-200 


.77 


.51 


.41 




-300 


.73 


.53 


.47 




+300 


.132 


.33 


.40 


400 


-100 


.76 


.20 


.14 




-200 


.72 


.47 


.39 




-300 


.685 


.51 


.46 




+300 


.166 


.35 


.40 


500 


-100 


.68 


.16 


.13 




-200 


.68 


.45 


.38 




-300 


.65 


.49 


.45 




+300 


.20 


.37 


.41 
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Figure Captions 



Fig. 1. The QCD cross-section for gluino pair production shown along with the electroweak 
cross-sections for chargino/neutralino production as functions of gluino mass for tan (5 = 
2 and /i = -100, -200, -300 and +300 GeV. 

Fig. 2. Electroweak cross-sections for chargino/neutralino production shown as functions of 
gluino mass for tan/3 = 10 and n = -100, -200, -300 and +300 GeV. 

Fig. 3. The gluino and the chargino/neutralino contributions to the LSD signal cross-section 
shown for the R- violating Yukawa couplings of eq. (5) at tan [3 = 2 and different values 
of /i. 

Fig. 4. The LSD signal cross-sections shown at tan/3 = 2 and /i = —100,-200,-300 and 
+300 GeV for various i?-violating Yukawa couplings. The solid, dot-dashed, long 
dashed and short dashed lines correspond to the couplings of eqs. (5), (6), (7) and (8) 
respectively. The trilepton signal cross-section for the last case is shown as crosses. 

Fig. 5. Same as Fig. 4 at tan/3 = 10. 
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